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Abstract : 

Long-period records of multiple Love waves from the 1957 
earthquake in Mongolia (Mg = 8.0) at Pasadena are analyzed and compared 
to synthetic seismograms, generated by the method of Kanamori (1970a, b). 

A fit in the time domain shows that the records are not consistent with 
the previous solution, achieved through a frequency-domain analysis of 
directivity by Ben-Menahem and Toksoz (1962). The solution asks for a 
shorter rupture of 270 km at a velocity of 3.5 km/s. The focal parameters 
are constrained by updating all the reported first-motion and are found 
to be : Strike = 103“, Dip = 53“ , Slip = 32“ . A seismic moment of 
1.8 10^® dynes-cm is obtained. These figures are also consistent with a 
time-domain analysis of Love waves at Palisades and Strasbourg, and of 
Rayleigh waves at Pasadena, with a directivity study of Love waves at 
Pasadena, and with static deformation and isosaismal data. A discussion 
is given of the relation between moment, magnitude and rupture area, and 
a comparison is made with other events in the same region: It is concluded 
that this earthquake does not exhibit an "intra-plate" behavior, but rather 
compares better with "inter-plate" events, such as the great Assam earth- 
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Introduction : 


At 03:37:45.0 TU, on December 4, 1957, a major earthquake 
(Mg = 8.0) occurred along the Bogdo fault in central Mongolia. The 
epicenter was subsequently computed at 45.3 ± 0.2 ®N by 99.3 ± 0.15 “E 
( Florensov and Solonenko , 1963). Although this event, known as the 
Gobi-Altal earthquake, did little civilian damage, due to its location 
under desert mountain ranges (Fig. 1), it was soon recognized as one of 
the major earthquakes in recent times, and an extensive field investiga- 
tion was conducted by a team of Russian and Mongolian scientists ( Floren- 
sov and Solonenko , 1963) . 

Furthermore, in view of the recent advent of plate tectonics, this 
earthquake could be regarded as the largest recent expression of intra-plate 
seismicity. It is therefore challenging to thoroughly investigate its 
mechanism. The only such attempt had been made by Ben-Menahem and Toksoz 
(1962) [hereafter referred to as BMT62], who studied multiple Rayleigh 
waves R^ to R^ recorded at Pasadena. However, at the time, the theory of 
surface-wave excitation was not fully developed and their study was focused 
on an investigation of directivity in the frequency domain. Their results 
asked for a fault rupture of 560 km, striking N100"E, at a rupture velocity 
of 3.5 km/s. The development of long-period surface-wave synthetic seis- 
mograms ( Kanamori , 1970a, b) allows a time-domain investigation of the 
rupture mechanism and warrants a complete reassessment of the focal and 
rupture process of the Gobi-Altai event. 

In Section 1, we discuss the focal solution of the earthquake, 
before determining the rupture process in section 2. A discussion of 
directivity effects is given in section 3. A computation of the seismic 
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moment Is given in section 4, and this figure Is discussed with respect 
to general Intra-plate tectonics. 

1. Focal process . 

The focal process of the Gobl-Altal earthquake has long been a 
challenge. At the time of the BMT62 paper, the solution could only be 
Inferred from surface deformation data. Hodgson et al .(1962) noted a high 
degree of inconsistency among the first-motion reports. They proposed a 
mechanism shown on Figure 2a, with strike (>!)), dip (6) and slip (X) 
angles of respectively 125.5, 55, and 24.3 degrees. Florensov and Solonen- 
ko (1963) show that their solution [ (j) = 107®, 6 = 41®, X = 40®] is rather 
poorly constrained, the azimuths of the planes, for example, being allowed 
to vary 35® around the chosen values (Figure 2b). They also report incon- 
sistencies in first-motion data. In view of the location of the station 
PAS (its azimuth from the direction of the fault varies from 60° to 95® 
depending on the model), it is crucial to achieve a better solution, if 
one is build synthetic seismograms. This is made possible by including 
mostly dilatational data from Indian stations (BOM, HYB, CCU, KOD), vhich 
were not used by Florensov and Solonenko, and mostly compressional data 
from Central Asia (DSB, KIR, AND, LAH) , which were not used by Hodgson 
et al . 

Inclusion of these data points yields a better-constrained 
solution: (j) = 103®, 6 = 53° and X = 32°, shown on Figure 2c. The rather 

large number of inconsistencies, mostly from European stations, may 
be due to the fact that most European stations lie in the direction oppo- 
site to faulting, and therefore give rise to emergent records. The 
long-period Galitsin record at DBN clearly shows a compression on Z, 
and a SW first motion on the horizontal instruments. Similar records at 
FIR also agree, as do both short- and long-period records at STR. Records 



at UCC show a self- inconsistent (Down, South, West) first motion. 

However, the shape of the signal is very similar to that at DBN (the 
two stations are only 150 km apart) and it is believed that the polari- 
ty of the Z instrument at UCC was wrong. This set of data clearly puts 
Western Europe into a compressional quadrant. The data from India and 
Central Asia constrains the solution to within a few degrees. It should be 
noted that the proposed solution, which represents a left-lateral strike- 
slip with a fair amount of thrust (see figure 2d), is consistent with 
the surface-rupture studies quoted from Florensov (1958) in BMT62. 

The solution is also basically consistent with the general pattern of 
Asian Tectonics, as described by Molnar and Tapponnier (1975). 

2. Rupture process from time-domain surface wave studies . 

Records used in this study were principally those obtained at 
’’AS, by the high-pass "30-90” 3-component seismometer ( Press et al . , 1958). 
and by the low-gain strainmeter. Additional Press-Ewing records from 
PAL and Galitsin records from STR were also used. Table 1 lists some 
parameters of the stations and records used in this study. The particular 
azimuthal location of PAS with respect to the fault elements is favorable 
to Love rather than Rayleigh waves. 

Love waves_at_Pasadena: Records of the three components of the 30-90 

instrument were digitised from 06:4A to 14:00 TU, and the horizontal 
ones were rotated into SH polarization. The phases G^, G^, and G^ were 
isolated and equalised to a station distance of 90 degrees. Components 
with periods lower than 60 sec. were at the same time removed. No attempt 
to process G^ and Gy was made, although these phases were readily ident- 
ifiable on the records (see figure 3). The computation of the synthetic 


0 

(' 


jeismograms followed the procedure described by Kanamorl and Clpar (1975). 
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C Several lengths of rupture were used, varying from 200 km to 560 km, 

this upper limit being the rupture length proposed by BMT62. A rupture 
velocity of 3.5 km/s, similar to that used in BMT62, was used throughout 
the study. A step-function type of source was chosen, and no attenuation 
of the source along the rupture was assumed. In section 3, we will discuss 
part of this last hypothesis on the basis of some directivity data. The 
azimuth of propagation of the faulting was taken to coincide with the 
trace of the fault pla le ( ((i = 103" ). 

Synthetics for G^, G^ and are shown on Figure 4, for fault 
lengths of 560 km (suggested by BMT62) and 270km (which turned out to give 
the best fit to our data), and are compared with the observed traces. It 
is readily seen that a fault length of 560 km cannot account for the 
observed wave shapes: The frequency content of the synthetic and G 5 
is clearly different from that of the observed phases, and the shape of 
the synthetic G^ exhibits a doublet, due to a rupture propagation backw- 
ards from the departing wave. On the other hand, a length of 270 km 
gives a most satisfactory fit to G 3 , and also to G 5 and G^, although the 
amplitude of G^ remains somewhat too large. 

Rayleigh waves at Pasadena. 

Synthetics were computed for the phase at PAS. It was found 
that the len^jth of the rupture has a rather limited Influence on its 
waveshape and that R^ cannot help in discriminating between various rup- 
ture lengths. The horizontal component of the phase R 2 , recorded on the 
low-gain strain Instrument, was also used: This system is equivalent to 
a mechanical seismometer with Tg = 90 sec. Figure 5 compares the observed 



7 . 


wave, filtered at. T = 50 sec., with synthetics for both 270 and 560 km. 

As can easily be seen, the shape of R2 at PAS is inconsistent with the 
560 km fault length. This is important, because Rayleigh waves were the 
only data used in BMT62. 

Other records. 

Figure 6 shows a similar comparison between observed and synthetic 

seismograms for at PAL. The other records obtained were those at DBN, 

FIR, UCC, and STR; these stations are located at an azimuth extremely 

unfavorable to Rayleigh waves, and, indeed, those are not easily observed 

on the records. The best record /as that of G2 at Strasbourg (Figure 7). 

/discrir.inate 

Unfortunately, it does not help between rupture lengths, because of sta- 
tion azimuth. This point will be discussed in the next section. 

As a partial conclusion to this section, it can be stated that 
most of the records of both G and R waves are incompatible with a rupture 
length of 560 km, whereas all agree with a shorter length of 270 km at 
3.5 km/s. 

3. Frequency-domain investigation . 

Our results conflict with those from BMT62, which were obtained 
in the frequency domain. In order to examine the cause of this discrepancy, 
we analysed our data in the frequency domain, and made a directivity plot 
similar to that of BMT62, for G3/G4 at PAS. The theoretical aspects of the 
problem of the directivity function of a moving source have been described 
by Ben-Menahem (1961) , who showed that a spectral component of frequency f , 
leaving the source at an angle 0 from the rupture direction, is attenuated 


by an amount : 



( 

( 

(1) D ( e ) 


where I is the rupture length, v the rupture velocity and c the surface- 
wave phase velocity at frequency f . For two successive multiple surface 
^c-'es, such as and G^, at the same station, equalised to an identical 
dis 'ance, the spectral ratio is then given by: 

D ( 0 ) 

( 2 ) 

D ( e + TT ) 

The experimental value of this ratio for G-j/G^ at PAS is shown on Figure 
8 and is compared with theoretical models for both fault lengths. Should 
the fault be 560 km long, a peak in the directivity function would be 
expected around 200 sec. , for which there is no evidence in the data. On 
the other hand, a reasonably good agreement is found with the theoretical 
spectrum ratio for 270 km. This indeed is a conclusion similar to that in 
BMT62, in which the authors noticed the disapearrance of the expected 
first Rayleigh wave directivity peak. Although the gain of the instrument 
is rather low at these frequencies, both the present investigation and 
Figures 6 & 7 in BMT62 show that the amplitude of the spectrum of R^, R 4 , 

G-j and G^ is indeed maximum around 200 sec. It is therefore not possible 
to account for the non-existence of the 200 sec. peaks through errors intr- 
oduced by the instruments. 
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Similarly, Figure 9 shows the spectrum of the R 2 record at PAS, 

on the long-period strain instrument. There is no evidence for the hole at 
/ required 

T = 195 sec. by the 560 km solution. The hole for 270 km falls in a rather 
rugged portion of the spectrum, for which no definite comment can be made. 

In the case of stations lying close to the azimuth opposite to the fault, 
such as most oi the European stations (see Table 1), formula (1) shows 
easily that even-order waves, leaving the source in the direction of rupture 
will be affected by directivity only in the higher frequency part of 
their spectrum. The only odd-order wave record of good quality from Europe 
was that of at STR. However, its Fourier spectrum is rather rugged and 
does not indicate any definite directivity effect. 

It is thus apparent that in the frequency domain, the analysis 
of directivity effects on Love and Rayleigh waves is not always a very 
powerful method of investigation of rupture mechanism. Difficulties in 
using this method were originally pointed out by Aki (1966). However, 
even in the Fourier space, most of the results from this study do 
agree better the shorter (270 km) length of the fault rather than with 
the longer (560 km) length suggested by BMT62 . We were not able to find 
any evidence for the fundamental hole expected around 200 sec. in G 2 ^ 
or R2n’ ^ km-length model. We would also like to point out that 
BMT62 used a data set of 2160 sec, which they extended to 20,000 sec. 
by adding zeroes. Therefore, their interval of sampling in the frequency 
domain cannot be shorter than 4.6 10 ** Hz = 1/(2160 sec.), and any 
feature of their spectra not extending over such an Interval has no real 
physical meaning. It shoud be noted that many of the spikes in Figure 11 
of BMT62 do not meet this criterion. Our analysis used a FFT over data 
samples of 2048 sec., and no substantial number of zeroes were added. 
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It therefore appears that the first peak of the directivity function, 
which Ben-Menahem and Toksoz identified as the //2-harmonic for a 560 kra 
rupture, is in fact the fundamental for a 270 km fault. 

As a conclusion to this section, it should be noted that the 
directivity plot on Figure 8 clearly shows strong spikes, extending over 
3 orders of magnitude. This restricts the rupture process and helps 
reject both a bilateral fault and an attenuating source function, as was 
discussed in BMT62. Also, the aftershock pattern does not suggest a bilateral 
faulting. However, there does remain a trade-off between velocity of rup- 
ture and length of fault, which can be resolved neither in frequency 
nor in time domain. The discussion of some static data will put some 
further constrains on the possible solution. 

4. Seismic moment and discussion : 

The comparison of amplitudes between synthetic and observed 
seismograms allows a computation of the seismic moment of the earthquake. 

records at PAS yield a value of 1.8 10^® dynes-cm; at PAS yields 

1.4 10^® dynes-cm; records at PAL ask for 2.3 10^® dynes-cm, although 

an accurate value of the gain of the Instrument at this time was not 
available; the G 2 record at STR yields 1.4 10^® dynes-cm. A reasonable 

value is therefore (1.8 ± 0.4)10^® dynes-cm. It is interesting to compare 
this figure with the static data reported in Florensov and Solonenko (1963). 
These authors report average static dislocations along the fault on the 
order of 8 m. The depth of the hypocenter is estimated by them at 16 km. 
Continuous surface breakage asks for a vertical extent of at least 16 km, 
but probably not greater than 40 kra. The width of the fault is then 20 to 
50 km, due to the Inclination of the fault plane, and the static moment 
M = liSD is found to be 1.3 to 3.2 times 10^® dynes-cm, for a fault 
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length of 270 km. In agreement with the long-period seismic data. 

?‘*gure 10 is adapted from Florensov and Solonenko (1963) 

[their figure 89, p. 251 of the English translation]. It shows the dis- 
tribution of the major aftershocks of the earthquake, and the line of 
isoseismlc degree 5. Superimposed on it, we have plotted the 270 km and the 
560 km faults. As can easily be seen, the 560 km fault would go well out 
of the isoseismlc line. Although the data about this line was obtained 
in desert land and no continuation of it is available into China, its 
significance is not to be underestimated, in view of the extensive field 
work reported by Florensov and Solonenko (1963). On the same figure, it 
is easy to see that the repartition of aftershocks is confined to a length 
of about 300 km along the fault. Again, this is an indication supporting 
the 270 km hypothesis. 

The magnitude of the Gobi-Altal earthquake has been computed 
by a number of agencies. The Bulletin of the Bureau International de 
Slsmologle (Strasbourg) (1958 summarizes the surface-wave magnitudes 
reported by various observatories. A mean-value of Mg = 8.04 ± 0.45 
results. A redetermination from 20-sec. surface waves at PAS yields 7.98. 

A value of 8.0 has therefore been adopted here. It then becomes possible 
to compare the parameters of this earthquake with those of others, in the 
way of the study by Kanamori and Anderson (1975). The preceding sections 
suggest the following values; Magnitude Mg = 8.0 ; moment M^ = 

1.8 10^® dynes-cm; rupture surface S = 10,000 km^ 

From these figures, it is possible to compute a stress drop 
Ao = 26 bars, and an apparent stress no = 10.5 bars. Although this last 
figure is not very accurate (and could probably vary from 3 to 30 bars), 
due to uncertainties in the magnitude-energy relation, it is still way 
under the average "intra-plate" apparent stress of 60 bars, determined 
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by Kanamorl and Anderson (1975). In this respect, the Gobi-Altai earth- 
quake plots among the inter-plate rather than Intra-plate family on 
Figure 4 of Kanamorl and Anderson's paper. Its next neighbors on the plot 
are the 1968 Tokachi-JkJ , 1969 Kurile, 1970 Peru, 1944 Tonankai and 1946 
Nankaido events, all major interplate earthquakes. The closest intraplate 
events are the 1952 Kern County and 1969 Portuguese earthquakes, which 
have much higher apparent stresses. 

A similar result may be obtained from the stress drop Ao. 

It should be emphasized that the relation between Mq and S, used in the 
determination of the stress drop, holds even for very large earthquakes, 
as pointed out by Kanamorl and Anderson. On Figure 2 of their paper, 
the Gobi-Altai earthquake is comparable to the 1952 and 1968 Tokachi-Oki, 
1966 Peru, 1969 Kurile, and 1944 Tonankai events. Again, intraplate events 
of comparable importance exhibit larger stress drops, and there seems to 
exist a major difference in terms of source parameters between the Gobi- 
Altai earthquake and other intra-plate events. In order to investigate this 
point furhter, we studied two other large events, which occurred in Mongolia 
and Sin-Klang ( Okal , 1976, to be published) : The Urumtchl earthquake of 
13 November 65, and the Khangal earthquake of 5 January 67 (see fijure 1). 
The parameters of these events are listed on Table 2. Purposedly, they 
were chosen so that the seismotectonics involved in them are different 
from that of the Gobi-Altai earthquake: The Urumtchi earthquake is a 
thrust event between the Tien-Shan fold belt and the Dzungarian block; 
the Khangai event is a pure strike-slip presumably along a North-South 
fault ( Moinar and Tapponnier , 1975), and possibly linked with the Lake 


Baikal rift system ( Florensov , 1968). We also include In Table 2 the 
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resjlt:s of Ben-Menahem et al .'s(197A) study of the great Assam earthquake 
of 15 August 50. Unfortunately, no data on fault dimensions is available 
for the Khangai and Urumtchi events, and only an apparent stress study can 
be made. Nevertheless, as can be seen, only the Khangai event clearly 
qualifies for intra-plate behavior. The Urumtchi earthquake has interme- 
diate properties, and the Gobi-Altai event joins the Assam earthquake into 
the inter-plate family. 

The interpretation of this feature is probably that the Bogdo 
fault, on which the Gobi-Altai event took place, is a majo’* expression of 
present Eurasian tectonism, along which a major left-lateral slip of China 
■ ^ -'ards the Pacific Ocean is taking place, as a result of the Himalayan 
collision, as proposed by Molnar et al . (1973) and Molnar and Tapponnier 
(1975), and possibly giving birth to a new "Chinese" plate. 

This indeed raises the question whether there can exist any major 
(Mg ^ ■•'.0) intra-plate event, or whether such an event is large enough as 
to make the two sides of its fault separate entities on a world-wide 
scale. A very challenging project in this respect would be the investi- 
gation of the huge 1905 Mongolian events ( Mg = 8.3 and 8.7 (?)) along 
the Bolnai fault. This might be done in the future, although the quality 
of the records involved is poor and a lot of uncertainty reigns over the 
characteristics of some of the instruments used at that time. 
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Cone jsion : 

We have shovm that a time-domain analysis of surface waves from 

the Gobi-Altai earthquake suggests a fault length of 270 km, which 

is much shorter than previously determined. Results in the frequency domain 

suggest the same. The focal parameters ( M = 1.8 10^® dynes-cm; M = 8.0) 

o s 

obtained from this study indicate that this earthquake has the characte- 
ristics of an inter-plate event, this being compatible with the present 
understanding of Chinese tectonics. 
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List of stations and records used in the study of the Gobi-Altai Earthquake. 
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Source parameters of Asian Earthquakes used in the present study. 
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Figure Captions. 


Figure 1. 


Figure 2. 


Figure 3. 


Figure 4. 


Figure 5. 


Figure 6 . 


Figure 7. 


Map of Mongolia, showing locations and focal mechanisms 
of the Gobl-Altal (4 Dec 57) earthquake and of the two 
events mentioned In Section 4. 

(a) Focal sphere plot of first-motion data obtained for 
the Gobl-Altal event by Hodgson et al . Black circles are 
compressions, open ones dilatations. 

(b) Focal mechanism achieved by Florensov and Solonenko. 
Heavy lines Indicate probable solution; dashed lines 
show possible uncertainty range. Black and open circles 
as In (a). 

(c) Focal mechanism obtained in this study. Black and open 
circles as in (a). 

(d) Scheme of focal dislocation. The arrow shows direction 
of movement of the hanging wall relative to the foot wall. 

Transverse (SH-polarized) trace of the Press-Ewing record 
at PAS. 

Observed and synthetic records of Love waves an PAS, 

(a) Synthetic seismograms for a 560 km rupture length 

(b) Synthetic seismograms for a 270 km rupture length 

(c) Observed signals. The 560km G 4 synthetic has been 
shifted upwards to enhance clarity. 

Observed and synthetic seismograms for the horizontal 
component of R 2 at PAS. (a),(b) and (c) as in figure 4. 

Observed and synthetic seismograms for G 3 at PAL. 
(a),(b),and (c) as in figure 4. 

Observed and synthetic seismograms for G 2 at STR. 

(a),(b) and (c) as in figure 4. 
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8. Spectral ratio G^/G^ at PAS. (a) Theoretical curve for 

a 560 km rupture; (b) Theoretical curve for a 270 km rupture; 
(c) Experimental curve. 

9. Spectral amplitude of observed R 2 trace at PAS. 

No hole in R 2 is observed at 195 sec., as would be 
required by the 560 km solution. 

10. Map of the epicentral area of the Gobi-Altai earthquake. 

The big star shows the epicenter; small dots indicate after- 
shocks; the rack line represents the isoseismic line of 
degree 5; the light trace is the 560km rupture fault sug- 
gested by BMT62 ; the heavier trace (displaced South to 
enhance clarity) is the 270 km fault propsed here. 
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